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Abstract. We measure feedback (heating rates) in damped Lya systems (DLAs) 
from the cooling rate of the neutral gas. Since cooling occurs through [C II] 158 
emission, we infer cooling from C II* 1335.7 absorption lines detected with 
HIRES on the Keck I telescope. The inferred heating rates are 30 times lower than 
for the Galaxy ISM. At z ~ 2.8 the implied star formation rate per unit area is 
logi/)*=— 2.4±0.3 Mq y~ x kpc~ 2 and the star formation rate per unit comoving 
volume logdp»/df=— O.8±O.2M0 y _1 Mpc -3 . This is the first measurement of star 
formation rates in objects that typify the protogalactic mass distribution. 



1 Introduction 

The purpose of this paper is to discuss a new method for probing physical 
conditions in DLAs. In principle, this technique can be used to probe their 
masses, sizes, and baryonic content. The focus of this paper is on star forma- 
tion. Specifically, in this paper I infer the star formation rates per unit area in 
DLAs by measuring the rate at which the neutral gas is heated. This is possi- 
ble because it is reasonable to attribute heating to far UV radiation emitted 
by massive stars. Since the cross-sectional area times the comoving density of 
the DLAs is known, we can derive the rate of star formation per unit comov- 
ing volume. Until now, comoving star formation rates have been obtained for 
highly luminous objects such as the Lyman Break Galaxies. Here I derive 
comoving star formation rates for objects that are more representative of the 
protogalactic mass distribution. 

2 [C II] 158 fim Emission From Damped Lya Systems 

[C II] 158 (im emission results from transitions between the 2 i"3/2 and 2 Pi/2 
fine-structure states in the ground term of C + . It dominates the cooling rate 
by the Galaxy ISM with a luminosity L([C II]) = 5xl0 7 L© f|. Most of 
the emission from the Galaxy and other nearby spirals arises in the diffuse 
neutral gas rather than from star-forming regions in spiral arms, or PDRs 
on the surfaces of molecular clouds (e.g. ||). The last point is especially 
relevant for DLAs where molecules are rarely detected ||. 

In order to estimate the [C II] 158 /im emission rates, Pottasch et al. [Q 
derived the following expression for the 158 fim luminosity per H atom: 

l c = N(CU*)his ul A ul /N(m) erg s" 1 H" 1 , (1) 
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where 7V(C II* ) is the column density of C + ions in the 2 P^/2 state, N(H I) 
is the H I column density, and A u i and hv u i are the Einstein A coefficient and 
energy of the 2 i"3/2 — > 2 Pi/2 transition. In fact, l c is the density— weighted 
average along the sightlinc of l cr (r) , the emission rate per H atom at a given 
displacement vector r. Here Z cr (r) = nc\i*{jc)A u ihu u i/riya{^)i where ncn* 
and tihi are the volume densities of H I and C II* . 

Figure 1 shows estimates of l c versus iV(H I), where 7V(C II* ) is deduced 
from the strength of C II* 1335.708 absorption and N(R I) from DLA ab- 
sorption. The small stars depict measurements for several sightlines through 
the ISM j|] and The large star depicts < l c >, i.e., l cr {r) averaged over 
the H I disk of the Galaxy. Figure 1 also shows our measurements of l c for 27 
DLAs. As neither 158 /jm nor 21 cm emission has been detected from any 
DLAs, < l c > cannot be inferred directly, where in this case < l c > corre- 
sponds to Z cr (r) integrated over the volume of the average DLA. Nevertheless, 
comparison between the two distributions of l c clearly indicates < l c > for 
the DLAs to be a factor of 30 or more lower than for the ISM. 
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Fig. 1. l c versus N(H I) for DLAs (solid circles) and sightlines through the ISM 
(stars). Upper limits are 2-cr. Large star is < l c >=L([C II])/(m_H M (H I)) where 
L([C II]) and M(H I) are the 158 luminosity and H I mass of the Galaxy. 
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3 PROPERTIES OF MULTI-PHASE MEDIA 
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Fig. 2. Keck HIRES absorption profiles for C IF 1335.7 and 3 resonance lines in 
the z = 3.0247 damped absorber toward Q0347— 348. Comparison between the C 
IF and unsaturated Fe II 1608 profiles shows evidence for a CNM at v — + 10 km 
s" 1 and WNM at - 10 km s" 1 

To determine the heating rate from measurements of l c we consider the 
thermal and ionization equilibria of gas subjected to heating and cooling. The 
result is a multi-phase medium in which the various phases are in pressure 
equilibrium. There arc indications that the H I gas in DLAs is in two phases: 
a cold (T > 50 K) neutral medium (CNM), and a warm (T - 8000K) neutral 
medium (WNM). The detection of 21 cm absorption in DLAs indicates the 
presence of a CNM since the 21 cm optical depth is inversely proportional to 
temperature Jj) . Evidence for a WNM comes from misalignment in velocity 
space between 21 cm absorption or C II* 1335.7 absorption on the one hand 
and UV resonance lines on the other (see Figure 2), since the CNM dominates 
158 /jm emission (Fig. 3c). 

To compute the properties of the multi-phase gas we focus on the thermal 
and ionization balance of the atomic CNM and WNM. By analogy with 
the ISM we assume that heating is dominated by photoelectric ejection of 
electrons from grains illuminated by far UV radiation (ft,zy=8-^13.6 eV). We 
adopt the formalism of |7j who considered small grains and PAHs and found 
the photoelectric heating rate at r to be given by 



r d (r) = 1.0xl(r 24 eG (r) ergs s" 1 IT 1 



(2) 
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In the last equation Gq(t), the radiation field incident at r, equals 47rJ(r) and 
is in units of Habing's [|| estimate of the local interstellar value (=1.6 x 10~ 3 
ergs cm~ 2 s^ 1 ). Here J is the mean intensity integrated between 8 and 13.6 
eV, and e is the far-UV heating efficiency. Following we include ionization 
and heating by cosmic rays and assume that the dust-to-gas ratio equals the 
metallicity, Z/Zq. As a result Id(r) is proportional to Z/Z Q . 

An example of the computed 2-phase diagrams is shown in Figure 3. Here 
we assume Go = 1.7 which corresponds to the far UV radiation field in the 
ISM (see H). We also let the cosmic ray ionization rate £ = 1.8xl0~ 17 
(see 0). For comparisons with the DLA data we let [Fe/H] (= logZ/Z Q ) 
= —1.5, the mean metallicity of the DLAs ( |J). Since regions of thermal 
stability occur where <9(logP)/<9(logn) > 0, a stable 2-phase medium, com- 
prising the WNM and CNM, can exist between the local pressure minimum 
and maximum. Figure 3c shows that l cr (r) asymptotically approaches a max- 
imum at log n > 1.5 cm~ 3 in the CNM where it dominates the cooling rate. 
Comparison with Figure 1 further indicates that this / cr (r) is much lower 
than the < l c > observed in the ISM, but falls within the range of l c spanned 
by the DLAs. 




Fig. 3. (a)Pressure vs. H density, n. (b) Electron fraction vs. n. (c) Photoelectric 
(PE) and cosmic-ray (CR) heating rates vs. n. Lya , [C II] (i.e., l cr (r)), grain 
recombination (Rec), and total (nA) cooling rates vs. n. (d) Temperature vs. n. 
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4 Star Formation Rates 




[Fe/H] 



Fig. 4. Computed l cr (r) vs. [Fe/H] for log</>„ = -2.8,-2.1, and -1.1 M^y^kpc' 2 . 
Data for DLAs show l c versus [Fe/H]. 



The previous discussion suggests that the far UV radiation field in high- 
z DLAs is similar to that in the Galaxy. Because the far UV intensity is 
proportional to the rate at which massive stars form flo|| , the star formation 
rates should be likewise similar. More precisely, the similarity is between the 
star formation rates per unit area, since solutions to the radiative transfer 
equation for plane parallel slabs of dust and stars show that 3 V is proportional 
to the luminosity per unit area, [i v |jTl|| . Here we compute J„ at the center 
of a uniform slab of stars and dust, and then integrate solutions for J v to J 
with |ll]] 's fit to the far UV spectrum of the ISM radiation field. Finally we 
combine J with eq. (2) to compute the heating rates. 

The curves in Figure 4 plot the resultant Z cr (r) versus [Fe/H]. The l cr {v) 
are computed in the CNM limit, l cr (r)=rd(r), and for r corresponding to the 
central points of uniform slabs. The curves are parameterized by 3 different 
values of ip* , the star formation rate per unit area. Comparison with the DLA 
data reveals significant variations of l c at fixed [Fc/H]. This indicates that ip* 
is distributed at constant metallicity. The data also show tentative evidence 
for the predicted decrease in l cr {r) with decreasing [Fe/H] at an apparent 
upper limit to the star formation rate of log-0*~ — 2.12 M Q y -1 kpc~ 2 . How- 
ever, we cannot rule out entirely a possibility currently under investigation; 
namely, that the DLA sightlines pass through gas containing only WNM [jl2| 
in which case the if)* in Figure 4 are lower limits. 
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Fig. 5. Star formation rates per unit comoving volume computed for (2 m = 1,Ha = 
0,/i=0.5 cosmology. Circles denote luminous galaxies |h| and point with error bars, 
DLAs. Horizontal error bar depicts redshift range of DLA sample. 



The star formation rate per unit comoving volume is given by 



dp*/dt 



< V* > (A/A p )xdN/dz 



(3) 



(1 + zfi-cdt/dz) 

where < ip* > is the average star formation rate per unit area, A and A p are 
the average cross-sectional area and cross-sectional area projected on the sky, 
and dN/dz is the number of DLAs per unit redshift interval. Assuming log< 
^» >= -2.4±0.3 M Q y- 1 kpc~ 2 and dN/dz=0.2A at z = 2.8 [|l3), the median 
redshift of our sample, we find that in the plane-parallel limit (i.e., A = 
2Ap), logdp*/dt =-0.8±0.2 MQy~ 1 Mpc~ 3 . This result, plotted in Figure 5, 
is compared with estimates derived from objects detected in emission such as 
the Lyman Break galaxies |L4j . It is encouraging that the DLA measurement 
is not widely different from measurements made with independent techniques. 
The similarity could mean that the DLAs and Lyman Break galaxies are the 
same entity. More likely, our estimates of ^* and dp*/dt suggest the DLAs 
to be characterized by low specific star formation rates and high comoving 
densities relative to the Lyman Break galaxies. 
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